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ABSTRACT 

We have obtained high-resolution Hi observations of N44, one of the 
largest Hll complexes in the Large Magellanic Cloud. The distribution and 
internal motions of the Hi gas show dynamic effects of fast stellar winds and 
supernova blasts. Numerous Hi holes are detected, with the most prominent 
two corresponding to the optically identihed superbubbles Shell 1 and Shell 
2. The Hi gas associated with Shell 1 shows an expansion pattern similar to 
that of the ionized gas shell, but the mass and kinetic energy of the Hi shell 
are 3-7 times those of the ionized gas shell. The total kinetic energy of the 
neutral and ionized gas of Shell 1 is still more than a factor of 5 lower than 
expected in a pressure-driven superbubble. It is possible that the central OB 
association was formed in a molecular cloud and a visible superbubble was 
not fully developed until the ambient molecular gas had been dissociated and 
cleared away. This hypothesis is supported by the existence of a molecular cloud 
toward N44 and the fact that the apparent dynamic age of the superbubble 
Shell 1 is much shorter than the age of its OB association LH47. Accelerated 
Hi gas is detected at the supernova remnant 0523—679. The mass and kinetic 
energy in the associated Hi gas are also much higher than those in the ionized 
gas of 0523—679. Studies of interstellar gas dynamics using ionized gas alone 
are clearly inadequate; neutral gas components must be included. 
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1. Introduction 

N44, cataloged by Henize (1956), is a luminous Hll complex in the Large Magellanic 
Cloud (LMC). It contains an assortment of compact Hll regions, filaments, and shells of all 
sizes, as well as three OB associations, LH47, 48, and 49 (Lucke & Hodge 1970). As shown 
in Figure 1, N44 is dominated by a prominent shell around LH47 in the central region (Chu 
& Mac Low 1990). In the surroundings are: diffuse Hll regions and hlaments to the east; 
compact Hll regions encompassing LH49 to the south; a large faint shell to the west; faint 
filaments to the north; a luminous, compact H II region around LH48 on the northeastern 
rim of the main shell; and numerous small, bright, single-star Hll regions in the outskirts 
of N44. 

Previous studies of ionized gas in N44 have yielded many interesting results. The main 
shell and the large shell to the west are designated as Shell 1 and Shell 2, respectively, 
and their expansion patterns have been studied by Meaburn & Laspias (1991). Shell 1 has 
been modeled as a superbubble using energy input implied by the observed massive stellar 
content; it is shown that the expansion velocity of Shell 1 is much higher than expected 
(Oey & Massey 1995). Diffuse X-ray emission has been detected in N44, indicating the 
existence of 10® K gas (Chu & Mac Low 1990; Wang & Helfand 1991). The X-ray emission 
within Shell 1 and Shell 2 has been suggested to be generated by supernova remnants 
(SNRs) shocking the shell walls; the outward extension of X-ray emission at the southern 
periphery of Shell 1 has been interpreted as a breakout; and the bright diffuse X-ray source 
at & — T northeast of Shell 1 has been identihed as a new SNR (Chu et al. 1993; Magnier 
et al. 1996). Following the naming convention of using truncated B1950 coordinates, this 
remnant will be named SNR 0523—679. 

Interactions between massive stars and the interstellar medium (ISM) have led to 
the derivation of a multi-phase ISM model (McKee & Ostriker 1977) and have been used 
to explain the disk-halo connection in a galaxy (Norman & Ikeuchi 1989). The physical 
structure of N44 vividly demonstrates the dynamical interaction between massive stars and 
the ISM. Thus, N44 provides an excellent laboratory to study in detail the distribution, 
physical conditions, and relationship of the different phases of the ISM in a star forming 
region. 

Previous studies of N44 have concentrated on only the warm and hot ionized gas. In 
order to carry out a comprehensive investigation of the multi-phase structure of N44, we 
have recently obtained high-resolution H l observations of N44 with the Australia Telescope 
Compact Array (ATCA) and high-dispersion, long-slit echelle spectra of N44 with the 
4m telescope at Cerro Tololo Inter-American Observatory. In this paper we report these 
observations (§2), describe the interaction between Hi and Hll gases (§3), analyze the Hi 
gas associated with Shell 1 and Shell 2 (§4), and study the acceleration of Hi gas by a SNR 
(§5). A discussion is given at the end (§6). 
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2. Observations 
2.1. HI Observations 

The Hi data of N44 are extracted from the Hi aperture synthesis survey of the LMC 
made with the Australia Telescope Compact Array (ATCA), which consists of 6 22-m 
antennas. The details of this survey are given by Kim et ah (1997) and preliminary results 
are presented by Kim and Staveley-Smith (1997). 

In summary, the observations were made in four conhgurations: 750A on 1995 February 
23 to March 11, 750B on 1996 January 27 to February 8, 750C on 1995 October 15 to 31, 
and 750D on 1994 October 26 to November 9. These conhgurations each have hve antennas 
with a maximum baseline of 750 m. The combined conhguration has 40 independent 
baselines ranging from 30 to 750 m, with a baseline increment of 15.3 m. The resultant 
angnlar resolntion is IfO for the Hi images presented in this paper. The largest angnlar 
strnctnre that the images are sensitive to is ~ 0?6. 

The observing band was centered at 1.419 GHz with a bandwidth of 4 MHz. The band 
was divided into 1024 channels; these channels were re-binned into 400 channels after an 
on-line application of Hanning smoothing and edge rejection. The hnal data cnbe nsed here 
has a velocity coverage of 190 to 387 km s“^, and a velocity resolntion of 1.65 km s“^. All 
HI and Ha velocities reported in this paper are heliocentric. 

The data cnbe of N44 was extracted from the fnll mosaic of the LMC, which consisted 
of 1344 pointing centers. Approximately fonr pointing centers were inclnded in the data 
cnbe of N44. The phase and amplitnde calibrators were PKS B0407—658 for some fields and 
PKS B0454—810 for the others. The primary ATCA calibrator, PKS B1934—638 (assnmed 
flnx density 14.9 Jy at 1.419 GHz), was observed at the start and end of each observing 
day. These observations served both for bandpass calibration and flnx density calibration. 

The data were edited, calibrated, and mosaicked in miriad. The pixel size is 20" and 
the size of the N44 cnbe is 23 x 23 x 120. Snpernniform weighting (Sramek and Schwab 
1989) was applied to the uv data with an additional Gaussian taper. The resnlting data 
were then Fonrier transformed and mosaicked in the image plane. The maximnm entropy 
method was used to deconvolve this cnbe. The final cnbe was constrncted by convolving the 
maximnm entropy model with a Gaussian of FWHM 1' x 1'. The hnal rms noise, measured 
in line-free parts of the hnal cnbe, is 9K. 


2.2. Echelle Observations 

To examine the kinematic properties of the ionized gas in N44, we obtained high- 
dispersion spectroscopic observations in the Ha-|-[N ll] lines with the echelle spectrograph 
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on the 4m telescope at Cerro Tololo Inter-American Observatory on 1995 January 18. 
Using a post-slit Ha filter and replacing the cross-disperser with a flat mirror, we were 
able to observe both Ha and [N ll]AA6548, 6583 lines over a long slit. The 79 lines mm“^ 
echelle grating and the red long focus camera were used. The data were recorded with a 
Tek 2048x2048 CCD. The 24 /rm pixel size corresponds to about 0.82 A (3.75 km s“^ at 
the Ha line) along the spectral direction and 0'.'267 along the spatial direction. The spectral 
coverage was limited by the hlter width to 125 A. The spatial coverage was limited by the 
slit length to ~ 4', with some vignetting in the outer l!5. The slit width was 250 pm, or 
1'.'66, and the instrumental FWHM was 16.1T0.8 km s“^. 

The wavelength calibration and geometric distortion correction were carried out with 
the use of Th-Ar lamp exposures in the beginning of the night. The observations were 
sufficiently deep that the geocoronal Ha component was detected in each frame. The 
geocoronal Ha component (at zero observed velocity) provided an accurate and convenient 
reference for velocity measurements. The observed velocities were then converted to 
heliocentric velocities. 

The echellograms are shown in Figure 1. Only Ha line is presented in the echellograms. 
The [N ll] lines are detected at a much lower S/N ratio than the Ha line, hence they are 
not presented here. 


2.3. Optical Imaging 

In order to compare the neutral gas components to the ionized gas in the N44 region, 
we obtained optical Ha emission-line images with a CCD camera mounted at the Newtonian 
focus of the Curtis Schmidt telescope at CTIO on 18 February 1994 UT. The detector was 
a front-illuminated Thomson 1024 x 1028 CCD with 19p pixels, giving a scale of T.'835 
pixeD^ and a held of view of 31f3. A narrow-bandpass Ha hlter (Ac = 6561 A, AA = 26 
A) was used to isolate the emission, and a red continuum hlter (Ac = 6840 A, AA = 95 A) 
was used to obtain images of the continuum background (images in other emission lines 
were also obtained, but those will be discussed in another paper). Multiple frames were 
obtained through each hlter, amounting to total integration times of 1200s in Ha and 600s 
in the continuum. The data were reduced with IRAF0, and multiple frames were shifted 
and combined to obtain the images shown in Figures 1 and 3b. 


®IRAF is distributed by the National Optical Astronomy Observatories (NOAO). 
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3. Interaction between H l and H ll Components 

The integrated Hi map of N44 (Figure 2) shows two clear depressions that correspond 
to the interiors of the optically dehned Shell 1 and Shell 2, and peaks that are adjacent 
to Ha features. The relationship between the Hi gas and the Hll gas is clearly complex. 
The ATCA channel maps of Hi (Figure 3a) kinematically resolve the distribution of Hi 
along the line of sight, and allow more precise identihcation of physical structures in the Hi 
distribution. Comparing these iso-velocity maps of Hi to the Ha image (Figure 3b), we 
may examine the relationship between the H l and H ll gas in greater detail and accuracy. 

From the comparisons of Hi and Ha images, we have found indications of: (1) Hi holes 
and expanding shells, (2) acceleration of Hi, (3) compression of Hi, and (4) the opening of 
a breakout in Hi. These phenomena are described below in more detail. 

Hi holes and expanding shells. The Ha image in Figure 1 reveals the prominent 
superbubble Shell 1 and the fainter superbubble Shell 2 in N44. Shell 1, around the OB 
association LH47, has a size of 70 pc x 50 pc; while Shell 2, with no obvious OB association 
inside, is 60 pc across. The Hi channel maps at Vhei = 292 - 302 km s“^ (Figure 3) show 
clear Hi holes at the central cavities of Shell 1 and Shell 2. Projected within the Hi holes 
are high-velocity Hi clumps (~15 pc diameter) at 5’^22™45®, —67°56'00" (J2000) at the 
velocity range of Vhei ~ 228 ~ 245 km s“^(see Figure 3). These high-velocity Hi clumps 
are most likely the brightest parts of the approaching hemisphere of an expanding H l shell 
associated with the optical shell. The receding hemispheres of these shells expand more 
slowly. The expansion of these H l shells is better visualized in the position-velocity plots 
(or L — V diagrams) in Figure 4. These shells will be analyzed in detail in §4. 

Besides the H l shells with obvious optical counterparts, we see two H l holes, centered 
at Vhei = 297 km s“^, without known expanding shells in the optical. Inspected closely, the 
HI hole at 5'^22™00®, —67°51'50" (J2000) is associated with sharp, curved optical hlaments 
that appear to delineate a shell structure, while the Hi hole at 5’^21™07®, —67°51'45" 
(J2000) appears to be associated with only faint, irregular optical hlaments. Interestingly, 
as shown in Figure 5, both of these Hi holes encompass regions of bright diffuse X-ray 
emission reported by Chu et al. (1993). The relationship between the X-ray-emitting 
ionized gas and the H l gas will be discussed further in §6. 

Acceleration of Hi. The Hi channel maps in Figure 3 show that the bulk of Hi gas 
is detected in the velocity range of 295 - 305 km s“^. High-velocity Hi gas with AH up to 
—80 km s“^ is detected in a portion of the SNR 0523—679. This high-velocity gas must have 
been accelerated by the SNR. High-velocity Hi gas {Vhei = 233 km s“^) is also detected in 
the region between the SNR and the superbubble Shell 1. This gas might consist of two 
components, one accelerated by the SNR and the other by Shell 1, as the Hi contours have 
an excellent correspondence with an Ha blister on the periphery of Shell 1 . 
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Compression of Hi. In the channel maps near the bulk velocity, Hi peaks are 
seen in regions between expanding Hi shells. It is possible that these peaks represent 
compressions caused by two shells expanding into each other. One example is the region 
at 5'^22“40®, —67°54'40" (J2000) between the SNR 0523—679 and the superbubble Shell 
1; the Hi column density, iV(Hl) = 3.9 x 10^^ cm“^, is a factor of 2 higher than those 
of adjacent regions along the rim of shell.|^ Another example is the region at 5*^23™06®, 
5^21™58®, —67°55'00" (J2000) bounded by Shell 1, Shell 2, and the Hi hole to the north. 
The Hi column density, A^(Hl) = 3.5 x 10^^ cm“^, is a factor of 2 higher than those of 
adjacent regions along the rim of shell. It is unlikely that the Hi peaks are caused by 
a low ionization in these directions, since the ionizing source of Shell 1 is located at the 
center of the shell. These H l peaks must be real enhanced column densities as a result of 
compression between expanding shells. 

Opening of a Breakont in H l. A breakout in Shell 1 toward the south has been 
suggested by Chu et ah (1993), based on an extension of diffuse X-ray emission toward the 
group of Hll regions associated with the OB association LH49. This explanation has been 
conhrmed by the presence of high-velocity ionized gas and a lower plasma temperature 
in the breakout region (Magnier et al. 1996). The Hi channel maps at Vhei of 300 - 305 
km s“^ show a clear extension of the Hi hole to the southeast, coinciding with the X-ray 
extension. This extension of the H l hole is likely the opening through which the breakout 
occurs. Near the southeast end of the X-ray extension, we see an Hi clump peaking at Vhei 
= 285 km s“^ with a column density A^(Hl) of 2.8 xlO^^ cm“^, a factor of 3 higher than 
the surrounding region. This Hi clump might have been impinged upon and compressed by 
the flow of hot, ionized gas in the breakout. 


4. The Superbubble Structure of Shell 1 and Shell 2 

Previous studies of the superbubble structure of N44 have used observations of only 
ionized gas (e.g., Chu et al. 1993; Magnier et al. 1996). It is found that the observed kinetic 
energy of the ionized Shell 1 of N44 is at least an order of magnitude lower than the value 
expected in Weaver et ah’s (1977) superbubble model (Magnier et al. 1996). Our new Hi 
observation of N44 reveals the existence of a neutral gas associated with Shell 1 and possibly 
Shell 2 as well. The kinetic energy of the neutral shell may alleviate the aforementioned 
discrepancy in energies. Below we derive the dynamic parameters separately for the ionized 
and neutral components of Shell 1. The energetics of Shell 2 is discussed near the end of 
this section. 


®The Hi velocity profile of this region has an extended red wing up to Vhei ~ 320 km s“^. However, the 
channel map centered at 310 km s“^ shows an unrelated feature extending from this region to the northeast. 
The compression is referred to the dense clump at the systemic bulk velocity. 
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The expansion of Shell 1 is not uniform. Its systemic velocity, determined from the Ha 
velocities at the rim of the shell, is 295 - 298 km s“^. Referenced to this systemic velocity, 
the receding hemisphere of the ionized gas shell has an expansion velocity of 30 km s“^, and 
the approaching side 45 km s“^. We will analyze the physical parameters of the approaching 
and receding hemispheres separately. These parameters are summarized in Table 2. 

The Ha echelle observation of the central region of Shell 1 shows that the receding 
component is three times as bright as the approaching component. Using photoelectrically 
calibrated PDS scans of the Curtis Schmidt plates of Kennicutt & Hodge (1986) and 
together with our CCD Ha image, we derive an emission measure of 214 cm“® pc for the 
receding component and 71 cm“® pc for the approaching component. To derive the rms 
shell density, we have assumed a uniform thin shell for Shell 1 to estimate its shell thickness. 
The Ha surface brightness prohle of Shell 1 shows a peak-to-center ratio of 13, indicating a 
fractional shell thickness AR/R of 0.02. The FWHM of the brightness peak at the shell rim 
is 0.1 times the shell radius. Given the nonuniform densities, the fractional shell thickness 
is most likely within the range of 0.02 - 0.1. We have derived the rms density, mass, and 
kinetic energy of the shell for shell thicknesses of 0.02 and 0.1, respectively, and listed them 
in Table 2. The rms density in the receding hemisphere is almost twice as high as that in 
the approaching hemisphere, but the kinetic energy of the receding hemisphere is 21% lower 
than that of the approaching hemisphere. The total mass and kinetic energy of the Hll 
shell are 2580 Mq and 3.5 x 10"^® ergs for a AR/R of 0.1, or 1210 Mq and 1.6 x 10^® ergs 
for a AR/R of 0.02. 

The HI material associated with Shell 1 of N44 shows a typical expansion structure in 
the position-velocity diagrams (Figure 4) along cuts crossing the central portion of Shell 1. 
Hi is detected in the velocity range of Vhei ~ 220 - 327 km s“^. The velocity gradient of 
the high-velocity component in the central cuts of Shell 1 indicates an expansion motion. 
At the center of Shell 1, the receding side of the Hi shell is expanding at a lower velocity, 
28±5 km s“^, than the approaching side, 50±5 km s“^. These expansion velocities are 
comparable to those of the ionized component of Shell 1. 

Following the analysis of the ionized shell, we discuss the receding and approaching 
hemispheres of the H l shell separately. The column densities, masses, and kinetic energies 
are listed in Table 2. The Hi column density of the receding hemisphere is 2.5 times that 
of the approaching hemisphere, but the kinetic energy of the receding hemisphere is 21% 
lower. This contrast is remarkably similar to that in the H ll shell. The total mass of the 
Hi shell is 8640 Mq, and the total kinetic energy in Hi is 1.1 x 10^° ergs. 

Note that the total mass and total kinetic energy of the H l component of Shell 1 are 3 - 
7 times higher than those of the Hll component. This implies that analyses of superbubble 
dynamics must include the neutral gas component! Note also that the asymmetric expansion 
of Shell 1, shown in both Hi and Hll components, indicates a stratihed interstellar medium 
with densities decreasing toward us. 



We may re-examine the energetics of Shell 1 with the additional H l information. The 
total stellar wind and supernova energy input implied by the observed stellar content (Oey 
& Massey 1995), the total kinetic energy derived in this paper, and the thermal energy of 
the hot interior derived from X-ray observations (Magnier et ah 1996) are summarized in 
Table 2. In a pressure-driven superbubble, the thermal and kinetic energies are expected 
to be 35/77 and 15/77 times the total mechanic energy input from stellar winds and 
supernovae (Weaver et ah 1977; Mac Low & McCray 1988). The expected thermal energy 
and shell kinetic energy are also listed in Table 2. The observed shell kinetic energy is 
dominated by that of the expanding Hi shell; still, there is at least a factor of 5 discrepancy 
between the observed and expected kinetic energies. The discrepancy might be caused by 
the breakout of Shell 1, since a large fraction of the stellar energy input may have been lost 
in the breakout. Other possibilities are discussed in §6.2. 

The superbubble Shell 2 of N44 appears as an H l hole in the channel maps as well 
as in the position-velocity diagrams (Figure 4). The receding side of the shell is detected 
along some cuts across the shell. The approaching side is much fainter; only a bright, 
high-velocity Hi knot is detected. (The rest of the approaching hemisphere is probably 
below our detection limit.) 

In a similar way, we describe the physical parameters of the Hi component of Shell 2. 
The receding part of the Hi shell has an expansion velocity of 30±5 km s“^, lower than that 
of the approaching side, 50±5 km s“^. The asymmetric expansion of the Hi in Shell 2 is 
similar to that of Shell 1; therefore. Shell 1 and Shell 2 probably share the same interstellar 
environment. The Hi mass of Shell 2 is 5850 Mq in the receding hemisphere and 3090 Mq 
in the approaching hemisphere. The total kinetic energy of the Hi in Shell 2 is 1.3 xlO^'^ 
ergs, 20% higher than that of the Hi in Shell 1. The stellar content of Shell 2 has not been 
as well studied as that of Shell 1, and the X-ray observations of Shell 2 have much lower 
S/N ratios, hence Shell 2 cannot be modeled in as much detail as Shell 1. 


5. The Physical Structure of the SNR 0523—679 

The SNR 0523—679, centered at 5'^23“06®, —67°53'15" (J2000), was hrst diagnosed by 
its diffuse X-ray emission (Chu et al. 1993). The diffuse X-ray emission is surrounded by 
curved Ho hlaments that are suggestive of a shell with a diameter of ~3'8. The Ho shell, 
called Shell 3 by Chu et al. (1993), is bounded by Hll regions on the west and south sides. 
The bright arc on the west side of the shell is also on the surface of the H ll region around 
the OB association LH48, indicating a physical interaction between the SNR and the H ll 
region. The southern part of the SNR shell is blended with a diffuse ring-like H ll region; 
however, it is not clear whether they physically interact. Within the boundary of the SNR 
shell, long Ha hlaments and dust lanes run from the NE corner to the SW corner, dissecting 
the shell into two lobes. The X-ray emission shows a corresponding 2-lobe structure, which 
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could be caused by the absorption in the dust lane. 

The internal motion of the SNR is revealed in the echellograms at slit positions S4 and 
S5 (Figure 1). The velocities at the shell rim sampled by these two slits are Vhei = 297 - 
304 km s“^, which will be adopted as the systemic velocity of the ionized SNR shell. The 
echellograms show very different velocity structures in the north lobe and the south lobe 
of the SNR. The north lobe shows prominent, blue-shifted material in the Ha line; the 
brighter parts have velocity offsets up to —125 km s“^, while the fainter parts have velocity 
offsets up to —150 km s“^. The receding side of the north lobe has much smaller velocity 
offsets, ~20 km s“^. The south lobe shows mostly red-shifted material in the Ha line, with 
velocity offsets up to -1-100 km s“^; however, blue-shifted material is also detected at the 
northern end of the south lobe, with velocity offsets up to —100 km s“^. 

For such an irregular pattern of expansion, it is impossible to determine the kinetic 
energy of the ionized gas in SNR 0523—679 using our limited echelle observations. We can 
nevertheless assume a uniform shell and derive upper limits on the mass and kinetic energy. 
The brightest part of the SNR along slit S5 has an emission measure of 25 cm“® pc. The 
apparent width of the SNR shell, 675 or 1.6 pc, provides an upper limit on the real shell 
thickness. The corresponding rms Ue is ~ 4 cm“^. The upper limits on the mass and kinetic 
energy are thus 370 Mq and 6 x 10^® ergs. 

The Hi gas associated with SNR 0523—679 has a very different velocity structure 
from that of the ionized gas described above. The bulk velocity of Hi gas toward the 
SNR is Vhei = 306±5 km s“^. Since this velocity is similar to the SNR’s systemic velocity 
determined with the Ha line, we will adopt this velocity as the systemic velocity of the H l 
gas associated with the SNR. The channel maps in Figure 3 show an Hi hole at Vhei = 297 
km s“^, although the size of the Hi hole is smaller than that of the Ha shell. 

The receding side of the H l gas associated with the SNR is clearly seen in the 
channel map at Vhei = 325 km s“^, suggesting an expansion velocity of 20-30 km s“^. The 
approaching side of the H l gas is more complex. Blue-shifted H l gas in the velocity range 
Vhei = 230 - 245 km s“^ is seen toward the SNR. However, this Hi gas also extends further 
south than the boundary of the SNR; the southern extension is particularly pronounced in 
the region between the SNR and Shell 1 at Vhei = 228 - 235 km s“^. As we explained in 
§3, Shell 1 might be responsible for the acceleration of the high-velocity Hi in the region 
between Shell 1 and the SNR. Interestingly, the highest-velocity iVhei = 226 km s“^) Hi 
component unambiguously associated with the SNR is located at the NE edge, instead of 
the center, of the optical shell. This might indicate that a breakout similar to that seen at 
the south edge of Shell 1 has occurred at this position. 

We may estimate the mass and kinetic energy of the H l gas associated with the SNR. 
The high-velocity Hi gas is more extended than the optical boundary of the SNR. To 
exclude the contribution from other sources, we use only the part of H l gas projected 



10 


within the optical boundary of the SNR. The mass of the approaching side of the H l gas is 
7460±1000 Mq and the kinetic energy of the SNR is 2.7 x 10^° ergs. The mass and kinetic 
energy of the Hi gas are much higher than those of the ionized gas in the SNR. This is 
similar to what we have found for the superbubble Shell 1. 

The thermal energy of the SNR can be derived from the ROSAT and ASCA 
observations in X-rays. The best Raymond and Smith (1977) model hts of these two sets 
of data give very different plasma temperature and absorption column density; however, 
the thermal energies derived from these two sets of parameters are similar. The best £t to 
the ROSAT Position Sensitive Proportional Counter data (Chu et al. 1993) gives a plasma 
temperature of kT = 0.835 keV, an absorption column density of log Nh = 20.68, and a 
normalization factor of log N‘^Vf /AnD"^ = 10.817, in the cgs units, where is the electron 
density, V is the volume, / is the volume filling factor, and D is the distance. The thermal 
energy of the SNR calculated from this model fit is 1.2 x 10^°(//0.5)^/^ ergs. The best 
fit to the ASCA Solid-state Imaging Spectrometers (SIS) data gives kT = 0.35 keV, log 
Nh = 21.48, and log N^Vf /= 11.326. The thermal energy from this model fit is 9.0 
X 10^®(//0.5)^/^ ergs. Considering the uncertainties of filling factor and spectral hts, the 
thermal energy of the SNR is probably (1.0±0.5)x 10^° ergs. 

The thermal energy in SNR 0523—679 is similar to that in each of the shells in the SNR 
DEML316 (Williams et al. 1997), despite the fact that SNR 0523—679 is in a much more 
inhomogeneous medium and has a much more nonuniform expansion. The kinetic energy in 
the ionized gas of SNR 0523—679 is of a similar order of magnitude as the thermal energy 
of SNR 0523—679. This behavior is similar to what is observed in DEM L 316. The total 
kinetic energy (Hi -|- ionized) of SNR 0523—679 is much larger than the thermal energy in 
the SNR’s hot interior. This behavior is similar to what is seen in N44’s superbubble Shell 
1 . 


6. Discussion 

6.1. Implications on the Hi Structure of the LMC 

The ATCA Hi survey of the LMC (Kim & Staveley-Smith 1997) shows Hi holes with 
sizes ranging from a few tens of pc to ~1400 pc. It is conceivable that these Hi holes 
have been created by the combined effects of fast stellar winds and supernova explosions. 
However, most of these holes do not have any recognizable stellar content or corresponding 
ionized gas shells, making follow-up studies difficult, if not impossible. N44, containing 
three OB associations at three different evolutionary stages, offers an ideal site for us to 
study how massive stars shape the interstellar medium. The current Hi structure in N44 is 
readily visible, whereas the ionized gas structure allows us to foresee the future Hi structure 
after the ionizing fluxes have terminated at the demise of all massive stars in N44. From 
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the study of N44, we hope to gain insight into the general Hi structure of the LMC. 

We hrst examine the validity of the apparent “holes” in the H l maps of N44. The 
integrated H l map of N44 shows two deep depressions. Both depressions have associated 
Hi shells, and both shells have optical counterparts, the superbubbles Shell 1 and Shell 
2. For Shell 1, the Hi shell and the ionized gas shell have similar expansion patterns and 
velocities: both have the far side receding at 30 km s“^ and the near side approaching at 
45-50 km s“^. For Shell 2, the Hi shell appears to expand faster than the ionized gas shell; 
the Hi shell has expansion velocities of 30-50 km s“^, while the ionized gas shell show only 
a velocity variation from 304 to 320 km s“^ (Meaburn & Laspias 1991). This different 
behavior might be caused by the different distribution of ionizing sources. Shell 1 has a 
central ionizing source, the OB association LH47, hence its neutral Hi shell is adjacent and 
exterior to its ionized gas shell. Shell 2, on the other hand, has no central OB association, 
and its ionizing flux might be provided by LH47’s O stars to the east of Shell 2. The 
structure of the ionized gas in Shell 2 thus depends of the relative locations of the ionizing 
stars, and would not have a one-to-one correspondence with the neutral H l gas shell. 

The channel maps of H l in N44 also show numerous holes. The most obvious H l holes 
are present in the map centered at the bulk velocity of N44, Vhei = 297 km s“^. In addition 
to the aforementioned superbubbles Shell 1 and Shell 2, two other holes are present to 
the north. As described in §3, one of these Hi holes has Ha hlaments delineating a shell 
structure, while the other has no evidence of a coherent shell structure in Ha images. 
These two H l holes must be real holes, as both regions show enhanced X-ray emission (see 
Figure 5), indicating that the holes are hlled with 10® K hot gas. It might be argued that 
the enhanced X-ray emission is resultant from a smaller foreground absorption. Since the 
integrated Hi map does not show obvious anti-correlation between Hi column density and 
X-ray surface brightness in N44, we consider this argument unlikely. These apparent H l 
holes in the channel maps must correspond to real cavities. 

Conversely, we next examine the H l maps for holes at regions where we expect them. 
The SNR 0523—679, having a diameter of ~ 4', is well resolved by our Hi maps. However, 
only a very shallow depression in H l column density is present in the integrated H l map or 
the channel map centered at 297 km s“^. This is perhaps understandable because the SNR 
0523—679 is in a low density medium, indicated by its fast expansion and large size. The 
HI gas that has been accelerated by the SNR also shows a patchy structure, consistent with 
a low density medium around the SNR. 

We may also expect H ll regions to show up as holes in H l maps. However, neither the 
H II region around LH48 nor the H ll region around LH49 shows a hole in the integrated 
or channel maps of H I. The lack of correspondence between young H ll regions and H l 
distribution indicates that the star formation must have taken place in molecular clouds 
and the Hll regions contain the dissociated and ionized molecular gas. This conclusion has 
been previous reached by Allen, Atherton, & Tilanus (1986) based on their Hi observations 
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of the spiral galaxy M83. The presence of molecular cloud in N44 is evidenced in the CO 
map by Cohen et ah (1988). CO maps made with the ESO-SEST at a higher resolution 
show molecular material associated with the B, C, and D components of N44, the compact 
H II regions along the west (the B component) and southwest (the C component) rim of 
Shell 1 and to the southeast (the D component, associated with LH 49) of Shell 1 (Israel 
et al. 1993; Chin et ah 1997). This molecular material must be the remnants of the natal 
clouds. 

Besides H l holes and expanding shells, we see clear evidence of acceleration of H l gas 
by superbubbles or SNRs, compression of H l gas between expanding shells, and breakout 
structures in N44. Similar phenomena must exist in other active star forming regions. It is 
conceivable that the interstellar medium has been shaped by multiple episodes of massive 
star formation, and that the complex Hi structure represents the cumulative effects of the 
previous massive star formation. 


6.2. Hi Shell and Superbubble Dynamics 

The most remarkable lesson we have learned from the H l study of superbubble in N44 
is that the neutral gas actually contains more mass and kinetic energy than the ionized gas! 
For example, in Shell 1 of N44, the mass and kinetic energy in the expanding Hi shell is 
3-7 times those of the ionized gas shell. Thus, the Hi gas is an important component in 
superbubbles. Previously, superbubble dynamics has been observed with optical emission 
lines, hence only the ionized gas shells have been considered. This clearly needs to be 
improved. 

It has been a standing problem that the observed kinetic energy in a wind-blown bubble 
is much too low compared to that expected from Weaver et al.’s (1977) pressure-driven 
bubble model (Oey & Massey 1995). We have summed the kinetic energy in the ionized gas 
shell and in the H l shell for the superbubble Shell 1 in N44, and found that the observed 
total kinetic energy is still at least a factor of 5 lower than expected. The observed ratio 
of kinetic energy to thermal energy is nevertheless closer to the theoretical value of 3/7, 
especially if the thermal energy in the breakout region is included. 

The superbubble dynamics of Shell 1 may have been signihcantly altered by the 
energy loss through the breakout; however, the breakout may not be the only cause for 
the discrepancy between the observed superbubble dynamics and theoretical predictions. 
The fact that the age of the OB association LH47, > 10 Myr (Oey & Massey 1995), is 
much larger than the dynamic age of Shell 1, < 1 Myr, indicates that the formation of a 
superbubble did not start as soon as the central OB association was formed. It is possible 
that LH47 was formed in a molecular cloud, and the currently visible superbubble Shell 
1 did not start to expand rapidly until it has broken out of the molecular cloud. This 



13 


hypothesis is supported by the remnant molecular clouds observed toward N44 (Israel et al. 
1993; Chin et ah 1997). Only the youngest Hll regions are still associated with molecular 
material, and these youngest H ll regions show neither H l holes nor H ll shell structures. 
Similar explanation has been proposed for the superbubble in Nil, and may be common 
among all superbubbles (Mac Low et al. 1997). 

Of course, off-center SNR shocks could have provided acceleration and contributes to 
the discrepancy between the observed superbubble dynamics and the theoretical predictions. 
Evidence of SNR shocks in N44 includes the excess X-ray emission (Chu et al. 1993) 
and strong [S II]/Ha ratios (Hunter 1994). Future observational studies of superbubble 
dynamics need to include both the neutral and ionized gas components and extend to 
superbubbles with a wide range of X-ray luminosities, so that additional SNR acceleration 
effects can be differentiated from normal superbubble dynamics. 
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FIGURE CAPTIONS 

Fig. 1.— Ha image of the N44 complex and echellograms. The echelle slit positions are 
marked in the Ha image. The slit ‘S3’ samples the center of the superbubble Shell 1. The 
slit “SI” samples the breakout region, or the X-ray South Bar (Magnier et ah 1996). ‘S4’ 
and ‘S5’ sample the lobes of the SNR 0523—679. The two narrow lines in the echellograms 
are the geocoronal Ha and the telluric OH line blend at 6577.284 A. 

Fig. 2.— An ATCA integrated HI map of the N44 complex. The optically defined ‘Shell F 
and ‘Shell 2’ are labeled on the HI contours. Corresponding HI holes in these positions can 
be clearly seen. 

Fig. 3.— The ATCA channel maps of N44’s HI datacube. The heliocentric velocity is given 
in the upper left corner. The H l contours are overplotted on (a) greay-scale images of H l 
themselves, and (b) an Ha image. The contour levels are 0.07, 0.10, 0.15, 0.20, 0.25, 0.30, 
0.35, 0.40, 0.45 Jy/Beam. 

Fig. 4.— The position-velocity plots of HI in the N44 complex. The contours corresponds 
to the associated pixel brightness > 2.5 a. The high-velocity gas associated with SNR might 
have unresolved gas component, which relate with an Ha blister on the periphery of ‘Shell 

r. 

Fig. 5.— Hi contours overplotted on an X-ray image of N44. The Hi contours are derived 
from the channel map centered at Vhei = 297 km s“^. The Hi contour levels are 0.07, 0.10, 
0.15, 0.20, 0.25, 0.30, 0.35, 0.40, 0.45 Jy/Beam. The X-ray image is a smoothed ROSAT 
PSPC image from Chu et ah (1993). 
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Table 1: Journal of Echelle Observations 


No. 

RA (J2000) 
(hh mm ss) 

Dec (J2000) 

( ° ) 

Slit 

Orientation 

Exposure 

(s) 

Remarks 

1 

05 22 32 

-68 00 06 

E-W 

1200 

South breakout region 

2 

05 22 32 

-68 00 46 

E-W 

1200 

40"S of slit 1 

3 

05 22 17 

-67 56 18 

E-W 

2x1200 

Center of Shell 1 

4 

05 23 01 

-67 52 02 

E-W 

2x1200 

SNR 

5 

05 23 05 

-67 54 09 

N-S 

2x1200 

SNR 

6 

05 22 32 

-68 00 26 

E-W 

1200 

20"S of slit 1 
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Table 2: The Physical Parameters of Shell 1 and Shell 2 


Parameter 


Shell 1“ 

Shell I'’ 

Shell 2 



Front 

Back 

Front 

Back 

Front 

Back 

The Hll hemisphere: 

Shell Thickness (pc) 


0.5 

0.5 

2.3 

2.3 



He (cm“^) 


11.9 

20.7 

5.6 

9.7 

— 

— 

Mass (M©) 


440 

770 

950 

1630 

— 

— 

Vexp (km s"^) 


45 

30 

45 

30 

— 

— 

Ekin (10^® ergs) 


0.89 

0.68 

1.9 

1.5 

— 

— 

The Hi hemisphere: 
N{HI) (10^0 cm-2) 


3.6 

1.4 



1.6 

3.0 

Mass (M©) 


6180 

2460 



3090 

5850 

Vexp (km s"^) 


28 

50 - 60 



50±5 

30±5 

Ekin (10^^ ergs) 


4.8 

6.1 



7.7 

5.2 

The whole shell: 

Hll Ekin (10"^® ergs) 


1.6 

- 3.4 





Hi Ekin (10^^ ergs) 



11 



13 


Total Ekin, observed (10^^ 

ergs) 

13 

- 14 





Eth, observed in shell (10^® ergs) 

3 

- 20 





Eth, observed in shell & breakout (10^® ergs) 

13 

- 30 





Wind and SN Energy Input (10^® ergs) 

400 

- 700 





Total Ekin, expected (10"^® 

ergs) 

78 ■ 

- 136 





Eth, expected (10"^® ergs) 


200 

- 300 






Notes: 

“ with an ionized shell thickness AR/R of 0.02 
^ with an ionized shell thickness AR/R of 0.1 









